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The pair correlation function obtained from the neutron diffraction data of Mozer, De Graaf, and 
Le Neindre, is given for liquid 4 He at several thermodynamic states above and below the superfluid 
transition. A method for smoothly and accurately extrapolating the pair function into the limit of zero 
internuclear separation is considered. The pair function is computed from an approximation integral 
equation for the ground state wave function and compared to the experimental results. 
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1 . Introduction 

The intensity of elastically scattered neutrons in an 
isotopic fluid defines [l] 1 the structure factor, S(/c;p, J), 
which is a function of k = \ k — k () | , where k {) is the wave 
vector of the neutron before scattering, and k the wave 
vector after scattering. The quantities p and T are 
respectively the density and temperature of the fluid 
sample. The structure factor is of particular interest 
as its Fourier transform, ^>{r\i\ p,T) , with respect to 
the internuclear distance r 12 , defines the local struc- 
ture of the fluid. 

g 2 (ri2;p,T) = l 

+ {27T 2 pr V y)~ l (Ikk sin /cr 12 [S(/c; p, T) — 1]. (1) 



The function #2(^12; p, T) is necessarily positive and 
is usually referred to as the pair correlation function; 
we will adopt this nomenclature. Both neutron and 
x-ray [2] measurements of S(k; p, T) in 4 He have been 
made, but data of high statistical accuracy, for several 
thermodynamic states in 1-3 K range, have only 
recently been reported [3|. The motivation for experi- 
ments at low temperatures is to explore the changes 
in the behavior of the local structure of liquid 4 He 
that might accompany the transition at T=T\ from 
normal to superfluid [4|. In addition, it is also known 
[5] that measurements of the isothermal density deriv- 
ative of £2(^12; P-* T) are a measure of the magnitude 
and range of the correlation between three atoms. 
Quite recently, Mozer, De Graaf, and Le Neindre 



[6, 7] have done neutron diffraction measurements of 
high statistical accuracy for several thermodynamic 
states in liquid 4 He. The thermodynamic states at 
which the scattering experiments were performed are 
listed in table 1. The table also includes the value of 
the superfluid transition temperature, 7\, for each 
density considered. All temperatures in table 1 are 
known to within ±0.005 K, and all densities to within 
±0.25 percent. 



Table 1. Thermodynamic states for g 2 (r; p, T) in liquid 4 He 

The bracketed number in the pressure column is the observed pressure in the Pascal unit; 
1 at m = 0.10 1325 M Pa. 



p (atoms/A 3 ) 


T(K) 


Tk(K) 


P(atm), 


(MPa) 


0.0230 


1.860 


2.135 


4.31 


(0.44) 


.0230 


2.050 


2.135 


3.96 


(0.40) 


.0230 


2.130 


2.135 


3.60 


(0.37) 


.0230 


2.140 


2.135 


3.60 


(0.37) 


.0230 


2.300 


2.135 


3.79 


(0.38) 


.0230 


2.840 


2.135 


5.25 


(0.54) 


.0236 


2.020 


2.107 


6.50 


(0.66) 


.0236 


2.860 


2.107 


7.81 


(0.79) 


.0245 


1.940 


2.050 


11.60 


(1.18) 


.0245 


2.860 


2.050 


12.51 


(1.27) 



1 Figures in brackets indicate the literature references at the end of this paper. 



This manuscript reports the values of #2(^12; P> T) 
obtained from applying eq (1) to the data of Mozer, De 
Graaf, and Le Neindre; we also present the results of 
two other computations involving the pair function. 
The first is an investigation of an algorithm for the 
short range (ri 2 — »0) behavior of #2(^2; p, T). The 
other computation gives a comparison of the experi- 
mental results for the pair function and the results 
obtained from an approximate integral equation for 
the ground state value, F=0, of the pair function. 
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Table 2. 



The pair correlation function, g (r ;p,T), for liquid ^He as 



r , temperature, T, and density, p. 



a function of internuclear distance, 



/°<r 3 > 


.0230 


.0230 


.0230 


.0230 


.0230 


.0230 


.0236 


.0236 


.0245 


.0245 


1.860 


2.050 


2.130 


2.140 


2.300 


2.840 


2.020 


2.860 


1.940 


2.860 


r, 2 (A) 






















.123 


.400 


.784 


1.714 


1.454 


1.762 


1.498 


1.242 


• 915 


.850 


.518 


.245 


.437 


.743 


1.559 


1.365 


1.612 


1.403 


1.184 


.928 


.826 


.574 


.368 


.481 


.681 


1.329 


1.228 


1.389 


1.254 


1.088 


.929 


.780 


.638 


.491 


.514 


.600 


1.060 


1.056 


1.126 


1.067 


• 961 


.899 


.706 


.674 


.614 


.518 


.509 


.789 


.865 


.857 


.857 


.810 


.823 


.606 


.658 


.736 


.485 


.415 


.549 


.670 


.612 


.645 


.643 


.700 


.484 


.579 


.859 


.417 


.324 


.361 


.485 


.407 


.448 


.475 


.543 


.355 


.449 


.982 


.325 


.242 


.228 


.320 


.248 


.280 


.318 


.372 


.233 


.294 


1.104 


.223 


.170 


.143 


.181 


.133 


.150 


.182 


.212 


.129 


.145 


1.227 


.127 


.111 


.089 


.071 


.051 


.059 


.075 


.080 


.049 


.025 


1.350 


.045 


.062 


.046 


-.011 


-.007 


.000 


-.001 


-.013 


-.007 


-.053 


1.473 


-.019 


.023 


.003 


-.068 


-.050 


-.035 


-.049 


-.069 


-.044 


-.090 


1.595 


-.064 


-.007 


-.045 


-.102 


-.083 


-.059 


-.075 


-.094 


-.066 


-.096 


1.718 


-.093 


-.027 


-.093 


-.118 


-.105 


-.079 


-.084 


-.099 


-.076 


-.082 


1.841 


-.108 


-.032 


-.128 


-.119 


-.115 


-.096 


-.079 


-.088 


-.072 


-.058 


1.963 


-.104 


-.020 


-.136 


-.105 


-.107 


-.101 


-.057 


-.060 


-.047 


-.026 


2.086 


-.073 


.015 


-.107 


-.071 


-.075 


-.083 


-.013 


-.012 


.004 


• 021 


2.209 


-.008 


.076 


-.033 


-.010 


-.011 


-.027 


.060 


.062 


.086 


.091 


2.332 


.100 


.167 


.084 


.087 


.090 


.076 


.166 


.168 


.200 


• 191 


2.454 


.252 


.289 


.241 


.226 


.233 


.228 


.307 


.307 


.344 


• 327 


2.577 


.442 


.440 


.428 


.404 


.412 


.420 


.476 


.473 


.513 


.494 


2.700 


.654 


.615 


.632 


.614 


.619 


.636 


.663 


.659 


.698 


.684 


2.823 


.867 


.803 


.838 


.834 


.835 


.853 


.853 


.851 


.887 


.882 


2.945 


1.060 


.990 


1.031 


1.043 


1.039 


1.053 


1.032 


1.035 


1.066 


1.071 


3.068 


1.216 


1.161 


1.197 


1.219 


1.213 


1.218 


1.187 


1.195 


1.223 


1.233 


3.191 


1.326 


1.299 


1.325 


1.345 


1.341 


1.339 


1.308 


1.320 


1.343 


1.355 


3.313 


1.390 


1.392 


1.40 8 


1.418 


1.418 


1.413 


1.390 


1.403 


1.420 


1.431 


3.436 


1.415 


1.435 


1.448 


1.444 


1.448 


1 .444 


1.430 


1.440 


1.451 


1.460 


3.559 


1.411 


1.430 


1.450 


1.435 


1.440 


1.438 


1.431 


1.436 


1.439 


1.447 


3.682 


1.387 


1.388 


1.421 


1.404 


1.407 


1.405 


1.398 


1.399 


1 .394 


1.403 


3.804 


1.347 


1.324 


1.371 


1.360 


1.358 


1.354 


1.340 


1.340 


1.327 


1.338 


3.927 


1.295 


1.251 


1.307 


1.307 


1.301 


1.293 


1.267 


1.269 


1.251 


1.263 


4.050 


1.230 


1.182 


1.237 


1.245 


1.236 


1.228 


1.190 


1.197 


1.175 


1.186 


4.172 


1.156 


1.122 


1.164 


1.174 


1.167 


1.161 


1.118 


1.129 


1.104 


1.112 


4.295 


1.078 


1.071 


1.092 


1.098 


1.095 


1.094 


1.056 


1.068 


1.041 


1.046 


4.418 


1.003 


1.028 


1.024 


1.021 


1.024 


1.029 


1.006 


1.014 


.988 


.987 


4.541 


.941 


.98 8 


.963 


.953 


.960 


.966 


.966 


.967 


.942 


.938 


4.663 


.898 


.949 


.911 


.900 


.908 


.910 


.932 


.925 


.904 


.897 


4.786 


.872 


.911 


.871 


.865 


.870 


.866 


.901 


.889 


.875 


.865 


4.909 


.861 


.879 


.843 


.846 


.846 


.837 


.874 


.862 


.853 


.842 


5.031 


.859 


.856 


.827 


.838 


.836 


.827 


.852 


.844 


.840 


.829 


5.154 


.860 


.845 


.823 


.836 


.834 


.831 


.840 


.836 


.836 


• 827 


5.277 


.861 


.846 


.830 


.837 


.839 


.847 


.839 


.840 


.842 


.834 


5.400 


.866 


.858 


.847 


.845 


.850 


.867 


.851 


.852 


.857 


.852 


5.522 


.876 


.879 


.871 


.860 


.867 


.886 


.874 


.873 


.880 


.878 


5.645 


.895 


.905 


.898 


.885 


.890 


.904 


.904 


.899 


.^11 


.910 


5.768 


.923 


.934 


.926 


.919 


.919 


.923 


.935 


.P29 


.944 


.944 


5.890 


.955 


.962 


.953 


.957 


.952 


.943 


.965 


.959 


.^78 


.978 


6.013 


.987 


.987 


.979 


.992 


.984 


.969 


.991 


.988 


1.007 


1.008 


6.136 


1.015 


1.009 


1.005 


1.021 


1.013 


.999 


1.013 


1.014 


1.030 


1.032 


6.259 


1.035 


1.027 


1 . 3.0 


1.040 


1.038 


1.030 


1.031 


1.035 


1.046 


1.050 


6.381 


1.049 


1.039 


1.056 


1.053 


1.057 


1.058 


1.046 


1.052 


1.057 


1.064 


6.504 


1.059 


1.048 


1.078 


1.061 


1.072 


1.078 


1.058 


1.065 


1.065 


1.074 


6.627 


1.067 


1.054 


1.094 


1.070 


1.082 


1.088 


1.068 


1.073 


1.071 


1.080 


6.750 


1.073 


1.059 


1.099 


1.079 


1.088 


1.089 


1.073 


1.077 


1.075 


1.083 


6.872 


1.076 


1.062 


1.094 


1.086 


1.087 


1.085 


1.074 


1.077 


1 .075 


1 .081 


6.995 


1.073 


1.064 


1.078 


1.087 


1.080 


1.076 


1.069 


1.073 


1.071 


1 .073 



We proceed next, to the results for g2(ri2; p, T) and 
then in section 3, we discuss the short range behavior. 
The computation of the pair function from the Percus- 
Yevick approximation will be given in section 4; the 
article concludes, at section 5, with a discussion of 
results. 



2. Pair Function from Scattering Data 

The results of applying eq (1) to the neutron diffrac- 
tion measurements are given in table 2, where we have 
listed £2(712; p, T) as a function of r i2 for all the 
thermodynamic states shown in table 1. 
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Table 2, continued. The pair correlation function, g 2 (r 12 ;p,T), for liquid ^He as 
distance, r , temperature, T, and density, p. 



a function of internuclear 



7.118 


1 .063 


1.062 


1.058 


1.079 


1.067 


1.066 


1.061 


1.065 


1.062 


1.061 


7.240 


1.047 


1.055 


1.038 


1.062 


1.049 


1.053 


1 .049 


1.053 


1.047 


1.045 


7.363 


1.028 


1.042 


1.022 


1.039 


1.030 


1.038 


1.034 


1 .037 


1.028 


1.027 


7.486 


1.009 


1.025 


1.011 


1.016 


1.012 


1.021 


1.01O 


1.020 


1.009 


1.008 


7.60O 


.095 


1.007 


1.003 


• 9Q5 


.998 


I.OO3 


1.003 


1.003 


.991 


.991 


7.^31 


.086 


• 991 


.096 


.982 


.988 


.987 


• °90 


.988 


.°77 


.977 


7.854 


.981 


.979 


.986 


.974 


.980 


.975 


.978 


.976 


.968 


.966 


7.977 


.979 


.972 


.074 


.970 


.973 


.^67 


.971 


.968 


.963 


.959 


8.099 


.975 


.970 


.963 


.968 


.966 


.963 


.966 


.063 


.961 


.956 


8.222 


.^71 


.971 


.955 


.965 


.960 


.962 


.064 


.°62 


.961 


.956 


8.345 


.967 


• .972 


.952 


.962 


.956 


.960 


.965 


.962 


.963 


.958 


8.46 8 


.965 


.973 


.955 


.959 


.955 


.959 


.967 


.963 


.965 


.963 


8.590 


.065 


.973 


.962 


.958 


.959 


.959 


.969 


.965 


.969 


.968 


8.713 


.970 


.973 


.970 


.961 


.966 


.961 


.972 


.969 


.974 


.974 


3.836 


.977 


.976 


.977 


.968 


.976 


.968 


.976 


.974 


.980 


• 980 


8.958 


.°85 


.981 


.983 


.977 


.985 


.978 


.980 


.981 


.987 


.987 


9.081 


.993 


.987 


.989 


.988 


.992 


• 991 


.087 


.988 


.994 


• QQ4 


9.204 


1.000 


.995 


.997 


.999 


.999 


1.003 


.994 


.997 


1.000 


1.002 


9.327 


1.005 


1.002 


1.006 


1.008 


1.005 


1.012 


1.002 


1.005 


1.007 


1.011 


9.449 


1.010 


1.007 


1.015 


1.016 


1.010 


1.018 


1.010 


1.011 


1.013 


1.018 


9.572 


1.C14 


1.011 


1.021 


1.022 


1.017 


1.020 


1.017 


1.017 


1.019 


1 .024 


9.695 


1.018 


1 .014 


1.025 


1.026 


1.023 


1.021 


1.022 


1.022 


1.023 


1.028 


9.817 


1.021 


1.016 


1.024 


1.027 


1.027 


1.022 


1.025 


1 .025 


1.026 


1.029 


9.940 


1.022 


1.019 


1.021 


1.026 


1.029 


1.022 


1.025 


1.026 


1.027 


1 .027 


10.063 


1.021 


1.020 


1.019 


1.024 


1.027 


1.022 


1.023 


1.026 


1.025 


1.023 


10.186 


1.018 


1.019 


1.018 


1.021 


1.022 


1.021 


1.019 


1.023 


1.021 


1.018 


10.308 


1.014 


1.015 


1.017 


1.017 


1.016 


1.018 


1.014 


1.018 


1.015 


1.012 


10.431 


1.009 


1.011 


1.017 


1.012 


1.010 


1.014 


1.009 


1.011 


1.008 


1.007 


10.554 


1.004 


1.006 


1.014 


1.007 


1.005 


1.008 


1 .004 


1 .005 


1.002 


1.001 


10.677 


1 .000 


1.002 


1.009 


1.002 


1.001 


1.003 


1.000 


.999 


.997 


.997 


10.799 


.997 


1.000 


1.001 


.998 


.999 


.998 


.998 


.996 


.994 


.994 


IO.022 


.^94 


.099 


.093 


.994 


.997 


.993 


.996 


.995 


.992 


.992 


11.045 


.992 


.999 


.987 


.991 


.994 


.989 


.905 


.994 


.091 


• POO 


11.167 


.991 


.998 


.985 


.990 


.991 


.986 


• °93 


.°93 


• °91 


.989 


11.290 


.990 


.095 


.985 


.989 


.987 


.985 


.991 


.991 


.990 


.987 


11.413 


.989 


.991 


.987 


.988 


.984 


.984 


.989 


.988 


.988 


.985 


11.536 


.989 


.988 


.988 


.986 


.982 


.985 


.988 


.985 


.987 


.984 


11.658 


.OO0 


.086 


.989 


.985 


.983 


• °88 


.087 


.084 


.086 


.984 


11 .781 


.992 


.987 


.988 


.985 


.986 


.992 


.988 


.085 


.987 


.987 


11.904 


.995 


.990 


.988 


.986 


.991 


• 995 


• 991 


.988 


.989 


.992 


12.026 


.998 


.995 


.990 


.990 


.996 


.99 8 


.994 


.993 


.994 


.997 


12.149 


1.000 


.999 


.993 


.994 


1.001 


1.000 


• °97 


.999 


.999 


1.002 


12.272 


1.002 


1.001 


.997 


.999 


1.003 


1.001 


1.001 


1.003 


1.003 


1.006 


12.395 


1.002 


1.002 


1.002 


1.003 


1.004 


1.002 


1.003 


1.005 


1.007 


1.007 


12.517 


1.002 


1.001 


1.005 


1.006 


1.003 


1.003 


1.004 


1.005 


1.008 


1.007 


12.640 


1.001 


1.001 


1.007 


1.007 


1.003 


1.005 


1.005 


1.005 


1.008 


1.006 


12.763 


1.001 


1.002 


1.008 


1.007 


1.003 


1.006 


1.005 


1.006 


1.007 


1.006 


12.885 


1.002 


1.003 


1.009 


1.007 


1.004 


1.006 


1.005 


1.007 


1.006 


1.007 


13.008 


1.004 


1.005 


1.009 


1.008 


1.006 


1.006 


1.005 


1.008 


1.006 


1.008 


13.131 


1.005 


1.006 


1.009 


1.009 


1.008 


1.005 


1.005 


1.009 


1.007 


1.009 


13.254 


1.006 


1.007 


1.008 


1.009 


1.009 


1.004 


1.005 


1.009 


1.007 


1.009 


13.376 


.1.005 


1.006 


1.006 


1.009 


1.009 


1.004 


1.005 


1.007 


1.006 


1.007 


13.499 


1.004 


1 .^04 


1.004 


1.007 


1.006 


1.003 


1.004 


1.005 


1.004 


1.004 


13.622 


1.002 


1.003 


1.002 


1.005 


1.004 


1.002 


1.003 


1.002 


1.002 


1.000 


13.744 


1.001 


1.002 


1.001 


1.002 • 


1.001 


1.001 


1.003 


.999 


1.000 


.998 


13.867 


1.000 


1.001 


1.000 


1.000 


.999 


1.000 


1.002 


.997 


.998 


.997 


13.990 


1.001 


1.001 


1.000 


.998 


.998 


.998 


1,001 


.996 


.997 


• 9Q6 


14.113 


1.001 


.1 .000 


.999 


.997 


.998 


.997 


.999 


.996 


.997 


.996 


14.235 


1.000 


.999 


.997 


.995 


.997 


.996 


.998 


.996 


.097 


.906 


14.358 


.998 


.997 


.995 


.994 


.997 


• °97 


• °96 


• °96 


.°97 


.996 


14.481 


.996 


.996 


.993 


.992 


.996 


.097 


.095 


.996 


.997 


.995 



There are two regions of n 2 which reflect the 
principal uncertainties in measurements of S(k; p, T); 
these are the small and large r V2 values. The structure 
factor is positive and must satisfy, 



and 



Km S(k; p, T) = pk B TK T 

^0 



lim S(k; p, T) = 1, 



(2) 



(3) 



where ks is the Boltzmann constant and k t the iso- 
thermal compressibility. The experimentally acces- 
sible regions of k, the magnitude of the change in 
the neutron wave vector upon scattering, is limited. 
Therefore, in practice, one must necessarily extrap- 
olate from a given small value of k to zero in order to 
satisfy eq (2) and one must also truncate S(k; p,T) at 
a finite value of k in satisfying eq (3). In the Fourier 
transform, the errors introduced in satisfying eq (2) 
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Table 2, continued. 



The pair correlation function, g 2 (r ; P> T), for liquid -He as a function of internuclear 
distance, r , temperature, T, and density, p. 



14.603 


.994 


.996 


.992 


.992 


.995 


.998 


.995 


.997 


.996 


.995 


14.726 


.995 


.996 


.993 


.994 


.994 


.998 


.995 


.997 


.996 


.995 


14.849 


.996 


.996 


.996 


.996 


.995 


.998 


.996 


.997 


.995 


.996 


14.972 


.999 


.998 


.999 


.999 


.997 


.998 


.997 


.997 


.996 


.997 


15.094 


1.001 


.998 


1.001 


1.002 


.999 


.999 


.998 


.997 


.996 


.998 


15.217 


1.002 


.999 


1.002 


1.003 


1.001 


.999 


.999 


.998 


.998 


.999 


15.340 


1.001 


.999 


1.001 


1.002 


1.002 


1.000 


.999 


.999 


1.000 


1.000 


15.463 


.999 


1.000 


1.000 


1.000 


1.002 


1.000 


1 .000 


1.001 


1.001 


1.002 


15.585 


.997 


1.000 


.999 


.999 


1.001 


1.000 


1.000 


1.002 


1.003 


1.003 


15.708 


.998 


1.000 


1.000 


.999 


1.000 


1.000 


1.001 


1.004 


1.004 


1.005 


15.831 


1.000 


1.000 


1.002 


1.000 


1.000 


1.000 


1.002 


1.004 


1.004 


1.005 


15.953 


1.002 


1.001 


1.003 


1.002 


1.000 


1.002 


1.002 


1.004 


1.003 


1.004 


16.076 


1.004 


1.002 


1.003 


1.003 


1.001 


1.004 


1.003 


1.004 


1.002 


1.003 


16.199 


1.005 


1.002 


1.002 


1.003 


1.001 


1.006 


1.003 


1.003 


1.002 


1.001 


16.322 


1.003 


1.003 


1.000 


1.003 


1.002 


1.006 


1.003 


1.002 


1.001 


1.000 


16.444 


1.001 


1.003 


.999 


1.001 


1.002 


1.005 


1.002 


1.000 


1.001 


1.001 


16.567 


.999 


1.002 


1.000 


1.001 


1.003 


1.002 


1.002 


1.000 


1.002 


1.001 


16.690 


.998 


1.001 


1.001 


1.001 


1.003 


.999 


1.001 


.999 


1.002 


1.002 


16.812 


1.000 


1.000 


1.003 


1.002 


1.002 


.996 


1.000 


.999 


1.001 


1.002 


16.935 


1.002 


1.000 


1.004 


1.002 


1.002 


.995 


1.000 


.999 


1.000 


1.001 


17.058 


1.003 


.999 


1.003 


1.002 


1.001 


.997 


.999 


1.000 


.999 


.998 


17.181 


1.002 


.999 


1.000 


1.001 


.999 


1.000 


.999 


.999 


.998 


.996 


17.303 


1.000 


.999 


.998 


.999 


.998 


1.002 


.999 


.999 


.997 


.995 


17.426 


.997 


.999 


.996 


.997 


.997 


1.003 


.999 


.998 


.997 


.995 


17.549 


.996 


.999 


.996 


.996 


.997 


1.001 


.999 


.998 


.997 


.997 


17.671 


.996 


.999 


.997 


.996 


.998 


.998 


.999 


.997 


.998 


.999 


17.794 


.998 


.999 


.999 


.998 


.999 


.996 


.999 


.997 


.998 


1.001 


17.917 


1.001 


.999 


1.001 


.999 


1.001 


.995 


.999 


.998 


.999 


1 .001 


18.040 


1.003 


.999 


1.002 


1.000 


1.001 


.996 


.999 


.999 


.999 


1 .001 


18.162 


1.003 


.999 


1.001 


1.000 


1.000 


.998 


.999 


1.001 


1.000 


1.000 


18.285 


1.001 


1.000 


1.000 


.999 


.999 


1.001 


.999 


1.002 


1.000 


.999 


18.408 


.998 


1.000 


.999 


.999 


.998 


1.002 


.999 


1.002 


1.001 
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effect the accuracy of the large ri 2 behavior of 
gi{t\2\ p,T). The errors introduced in satisfying eq 
(3) give the unphysical oscillations about the value 
zero which are shown in table 2. Mozer et al. have 
given [6, 7] an extended discussion of the extrapolation 
to k=0, and we will, in the following section, discuss 
an algorithm to smoothly and accurately extend 
g2 (ri2; p, 7") to the value zero which it must necessarily 
have at r 12 = 0. It suffices to mention here that in the 
numerical evaluation of the Fourier transform, eq (1), 
S(k;p,T) was truncated to unity for k^I A 1 . The 
statistical accuracy of the resulting pair functions is 
estimated to be within ± 3 percent. 



3. Short Range Behavior 

The negative values of g 2 {r i2 ;p,T) in table 2 are 
unphysical and are introduced into eq (1) by truncating 
S(k;p,T). For physical systems, the pair function 
must approach zero in value as the internuclear 
distance approaches zero, and the approach to zero 
is presumably monotone. 

We consider a method for smoothly extrapolating the 
results in table 2 to the value g 2 = 0; as we will discuss, 
the method also applies to classical liquids. Consider 
the pair function for r i2 ^ d and denote the function on 
this interval by g 2 (r i2 ^ d). The form [8, 9] 

g 2 (r 12 ^ d)=g 2 (d; p, T) exp (x[l - (^J"]} (4) 



will be investigated, where we choose d as the smallest 
value of r 12 for which the pair function has a statistical 
accuracy equivalent to that at larger r l2 values. The 
number g 2 (d; p, T) is the value of g 2 (r l2 ; p, T) at 
r i2 = d. For a given value of m, the quantity k is 
determined by the continuity of the derivative of the 
pair function and g 2 (ri 2 ^ d) , 



k = 



g'Ar^P, T) \ 
m \g 2 (r l2 ; p, T)/ ri2 = { 



(5) 



The values m = 5 and m = 10 were considered and 
the typical agreement is shown in figures 1 and 2. In 
figure 1, the solid curve is the experimental result for 
4 He at p = 0.0230/ A 3 and T= 2.050 K; the crisscross 
denotes the value of d. The dashed curve is for m = 5 
and k= 2.423, and the dashed-dotted curve for m= 10 
and k = 1.211. Similarly in figure 2, the solid curve is 
the experimental result for 4 He at p = 0.0230/ A 3 and 
T = 2.840 K; the dashed curve is eq (2) with m = 5, 
k = 2.032 and the dashed-dotted curve corresponds to 
fli = 10, k = 1.016. As we have already mentioned, the 
results in figures 1 and 2 are typical of the results 
obtained using all the helium data in table 2, and from 
this we infer m = 5 is a good choice for 4 He. To deter- 
mine whether the algorithm also has utility in classical 
liquids, we investigated the data of Yarnell, Katz, 
Wenzel, and Koenig [10] taken for liquid argon near 
its triple point, p = 0.0215/ A 3 and T = 85 K. The results 
are shown in figure 3 where the solid curve is the exper- 
imental data, the dashed-dotted curve is eq (2) with 
m = 12, k = 2.670 and the dashed curve corresponds 
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Figure 1. The short range behavior of g 2 (ri 2 ) in liquid 4 He. 

The solid curve is the experimental result at p= 0.0230/A 3 , T= 2.050 K; the dashed curve 
is g( n 2 =s d) , eq (4), with m = 5, A = 2.423; the dashed-dotted curve is the same function 
with m = 10, \ = 1.21 1. The crisscross denotes the value of d. 
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FIGURE 2. The same as figure 1, except the experimental result is 
at 1 = 2.840 K; the dashed curve is for m=5, X = 2.032; the dashed- 
dotted curve for m= 10, X= 1.016. 

For 2.332 A «£ r, 2 < 2.867. the dashed curve is identical to the solid curve. 



to m = 5, A = 6.407. Relative to the data for 4 He, the Shiff and Verlet [14] have found from computer 

argon results suggest that for classical fluids, m— 12 is simulation studies, that for a Hamiltonian whose in- 
a better choice. teraction energy, U(r 12) , is the Lennard-Jones function, 



4. Percus-Yevick Approximation 

For a fluid system at 7 1 =0, where the wave function, 
^°\ f° r a set of N atoms is given by the Bijl-Dingle- 
Jastrow form [11], 



*<?>= II exp{-i*(ry)}, 



(6) 



the Percus-Yevick (P.Y.) approximation [12], may be 
used [13] to obtain an approximate integral equation 
for g2 {r\2 ) in terms of (f> (rn ) . The result is, 

g2(ri 2 ) = exp{-(/>(ri2)}[l + p I ^1^2(^3) [gifa) — 1] 

[l-exp{<Mn 3 )}]]. (7) 



U(r 12 ) ~ 



ri2 



r-fer 



the function 



0(r 12 ) 



= (^) S ' 



(8) 



minimizes the ground state energy. Using eq (8) and 
the algorithm given by Mandel, Bearman, and Bearman 
[15], we have solved eq (7) for £2(^12). All computations 
were done in units of o - , which for 4 He is [14], 



0-= 2.556 A. 



(9) 
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FIGURE 3. The short range behavior of g2(ri 2 ) in liquid At. 

The solid curve is the. experimental result from reference [8J; p = 0.0215/A 3 , T= 85 K. 
The dashed curve is eq (4) with m = S, \ = 6.407; the dashed-dotted curve is the same 
function with m=12, X = 2.700. The crisscross denotes the value of d. For the largest r i2 
values shown, the dashed-dotted curve is identical to the solid curve. 
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FIGURE 4. The pair correlation function for liquid 4 He in units of 
o- = 2.556 A at per 3 = 0.409. 

The solid curve is the experimental result at T— 1.940K. The dashed curve is the P.Y. 
approximation, eq (7), with the ground state function defined by eqs (6) and (8). 



The results for pcr 3 = d 409, which with eq (9) cor- 
responds to p = 0.0245/ A 3 , are shown in figure 4. The 
dashed curve is the P.Y. solution and the solid curve 
is the measured value atp = 0.0245/A 3 and r= 1.940 K 
In considering the relative agreement, which is within 
± 15 percent, we wish to emphasize that there are three 
sources of error: (a) the uncertainty between eqs (6) 
and (8) and the true ground state of helium, (b) the 
P.Y. approximation, (c) the contribution of excited 
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FIGURE 5. The difference, gf xp '— gV'-, at pcr 3 =0.394. 

The experimental result is at 7=2.020 K. 
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FIGURE 6. Same as figure 5, except pcr 3 = 0.394 and experimental 
result is at T— 1.860 K. 




2.10 



FIGURE 7. Same as figure 5, except pa 3 = 0.365 and experimental 
result (ref. [14]) is at 1 = 0.79 K. 
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quantum states. Keeping this in mind, it is of interest 
to consider effects of finite temperature. The results 
arc- shown in figures 5, 6, and 7, where we have plotted 
g expt._ g p.Y. for po -3 = o.394, 0.384, and 0.365 respec- 
tively. The experimental data used in constructing 
figures 5 and 6 correspond, respectively, to the measure- 
ments at p = 0.0236/A 3 , T= 2.020 K and p = 0.0230/A 3 , 
T= 1.860 K. The experimental data used in figure 7 
were obtained from the x-ray measurements of Achter 
and Meyer [3]; these data are at p = 0.0219/A 3 (pa 3 
= 0.365) and T=0.19 K. The results in figures 5, 6, 
and 7 indicate that at intermediate values of the 
internuclear separation, the P.Y. approximation 
characteristically underestimates the pair correlation 
in liquid 4 He, and at larger values of ri 2 /<x, the P.Y. 
approximation characteristically overestimates the 
measured values. Little, if anything, can be inferred 
from the differences at small internuclear separations, 
fWcr <1, since the measurements are, as we have 
already remarked, not as accurate as the data at 
larger separations. 

5. Discussion 

In a separate article [16], we have employed the data 
in table 2 in the study of several properties of the pair 
and triplet correlation functions. We reiterate one 
property that is particularly germane to this article. 
We consider table 2 in the vicinity of the first maximum 
of £2(^12; p, T), r i2 ~ 3.5 A, and the second minimum, 
Hz ~ 5.3 A. The magnitude of #2(^12; p, T) — \ is re- 
ferred to as the structure in the pair function. Table 
2 shows that below 7\, the structure is a decreasing 
function of temperature as the temperature is lowered 
at fixed density. The same dependence on temperature 
is observed as the temperature is lowered, at fixed 
density, from above 7\, but below the boiling point, 
to temperatures below the superfluid transition. For 
classical fluids the pair correlation function is [16] an 
increasing function of temperature as the temperature 
is lowered at fixed density. 

The algorithm given in eq (4) is, from figures 1, 2, 
and 3, a simple and accurate procedure for smoothly 
extending the values of the pair function inside the 
first maximum, to the value gt — 0. The short range 
behavior is very difficult to measure and perhaps the 
algorithm will prove useful in presenting experimental 



results from diffraction experiments. As we have al- 
ready noted, the method appears to have utility in 
classical as well as quantum fluids. 

It should be emphasized that the computations given 
in section 4 strictly apply to a fluid whose wave func- 
tion at the absolute zero of temperature is given by eq 
(6). With this in mind, the results in figures 4 thru 7 
show the P.Y. approximation is with± 15 percent of the 
experimental results. The P.Y. approximation typically 
underestimates the magnitude and extent of the first 
peak in #2(^12) but overestimates the function at 
somewhat larger (r 12 "> 1.8a*) values of the internuclear 
separation. As expected, the agreement between the 
calculated and experimental results improves when 
the calculated results are compared to data at lower 
temperatures. 



The authors express their appreciation to B. Mozer 
for the use of the diffraction data prior to publication. 
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